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The issue of the parameterization of small scale ( "sub-grid" ) turbulence is addressed in the 
context of passive scalar transport. We focus on the Kraichnan advection model which 
lends itself to the analytical investigation of the closure problem. We derive systematically 
the dynamical equations which rule the evolution of the coarse-grained scalar field. At 
the lowest-order approximation in l/r, I being the characteristic scale of the filter defining 
the coarse-grained scalar field and r the inertial range separation, we recover the classical 
eddy-diffusivity parameterization of small scales. At the next-leading order a dynamical 
closure is obtained. The latter outperforms the classical model and is therefore a natural 
candidate for subgrid modelling of scalar transport in generic turbulent flows. 



1. Introduction 

One of the most striking characteristics of hydrodynamic turbulence is the presence of 
a wide range of active length and time scales. These scales are strongly and nonlinear ly 
coupled, a fact that makes analytical approaches, at best, impractical. The situation does 
not look better for direct numerical simulations of turbulent systems: to fully resolve a 
turbulent flow requires approximately (L/77) 3 / 4 grid points in each spatial direction (see, 
for example, Frisch 1995), L and rj being the integral scale and the dissipation scale 
respectively. In the atmosphere, for instance, the ratio L/rj may become of the order of 
10 10 (77 - 1CT 3 m and L ~ 10 7 m) thus requiring the dynamical description of 10 
degrees of freedom. This remains, up to now and probably also in the near feature, a 
prohibitive task. 

To overcome the problem, "coarse-grained" versions of the original hydrodynamic equa- 
tions are often considered in order to describe large-scale features of the original full sys- 
tem. The Large-Eddy Simulation (LES) technique is probably the most popular example 
(Meneveau & Katz 2000). The success of such strategy is however strongly dependent 
on the realism of the description of small scales in terms of the large, explicitly resolved, 
scales. The problem of representing small unresolved scales in the absence of scale sep- 
aration - the long-known closure problem - attracts a great deal of attention in many 
domains ranging from geophysics to engineering (McComb 1992), and is one among the 
many challenges of turbulence theory. 

Our goal here is to shed some light on this aspect within the context of scalar turbulence 
where considerable progresses have been achieved in the last few years (Shraiman & Siggia 
2000, Falkovich et al. 2001). For this purpose we will consider a particular model of scalar 
transport (Kraichnan 1968, Kraichnan 1994) where the LES strategy can be formulated 
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and the problem of relating unresolved scales to resolved ones can be successfully attacked 
analytically. 

In this respect, the Kraichnan model has some characteristics of paramount impor- 
tance: 

- Exact expressions for relevant statistical observables can be derived from first prin- 
ciples, that is from Eq. I|2.1[) : this amounts to saying that the observables for the "fully 
resolved case" are known. An example is the expression (|3.3[) for the second order scalar 
structure function, an observable tightly related to the Fourier spectrum of the scalar 
field; 

- Closures for the large-scale dynamics can be derived in a systematical way (sec §|5} , 
and their predictions can be analytically checked against the exact solution. 

Those features make the Kraichnan model an ideal playground for studying LES closures. 



2. LES closures for passive scalar turbulence 

Scalar transport is governed by the advection-diffusion equation 

d t + v -V6> = n Q A9 + f , (2.1) 

describing the evolution of a passive scalar field 9(x, t) - e.g. temperature when buoyancy 
effects are negligible - advected by an incompressible velocity field v(x,t). Scalar fluctu- 
ations are injected into the system at the large scale L by the forcing term /. Dissipation 
takes on at small scales 77 due to the molecular diffusivity kq. 

The coarse-grained scalar and velocity fields, denoted by 9 and v are obtained by 
convolving the original, fully resolved, fields with a filter Gi with characteristic scale I 
(77 < I < L) 

6(x,t) = J Gi(x-x')e(x',t)dx' , (2.2) 

v(x,t) = J Gi(x- x')v(x',t)dx' . (2.3) 
The equation for 9 derived from Eq. I|2.1|l is: 

d t 6 + v- V 9 = K A6 + f- (C + S\ (2.4) 
where C is analogous to the Leonard stress 

£ = v-\76-v-\76, (2.5) 

and S is defined as 

s = v 1 ■ ve + v' ■ V0' + v ■ ve' . (2.6) 

The small-scale fields 9' and v' are defined as 6' = 9 — 9 and v' = v — v. 

The purpose of LES closures is to express C and S in terms of the large scale fields 
9 and v. Once this goal is accomplished, Eq. I|2.4|l can be numerically integrated on a 
mesh of spacing I, rather than 77 as it should be required for the integration of the full 
system (|2.1ll . with an enormous gain in memory and CPU time requirements. 

Unfortunately, no general closed expression for C and S in terms of 9 and v is available. 
A remarkable exception is the case where there is a marked scale separation between 
velocity and scalar length and timescales. It is then possible to show (see, e.g. , Biferale 
et at 1995, Mazzino 1997) that the effect of unresolved scales is just the renormalization 
of the molecular diffusion coefficient kq to an enhanced eddy-diffusivity k°" (generally 
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speaking, an eddy-diffusivity tensor nf- ). General expressions of the eddy-diffusivity as 
a function of the flow properties do not exist, and in most cases n° a can be determined 
only numerically. 

Here, our aim is to consider the challenging situation where there is no scale separa- 
tion between velocity and scalar and explore, in such context, the existence of effective 
equations for 9. 

Our procedure to derive the closed coarse-grained dynamical equations is the following: 

(i) starting form first principles, that is from Eq. I|2.1[) . we take advantage of the 
peculiarity of the Kraichnan model to derive an exact, yet unclosed, statistical equation 
for the correlation function of the filtered scalar field (6(x,t)6(x + r,t)); 

(ii) the statistical averages appearing in this equation, which involve small-scale fields, 
are then expressed in terms of correlations of large-scale fields only, at a given order of 
approximation in l/r; 

(in) from the statistically closed equations we consistently infer the dynamical closed 
equations for 9. 

No supplementary assumptions are made in performing this procedure. Although the 
method relies heavily on peculiar characteristics of the Kraichnan model, we believe that 
the results are relevant to generic passive scalar turbulence as well. Our claim is supported 
by several numerical and analytical evidences gathered in the last few years showing that 
most of the phenomenology of scalar turbulence is captured by the Kraichnan model. 
For an exhaustive review on this aspect, see, e.g. Falkovich et al. (2001). 

To illustrate the power of our approach we anticipate here the main results of this 
paper, postponing their derivation to the following sections. 

Carrying out the procedure at the lowest significant order in l/r yields the following 
effective equation for the coarse-grained field: 

d t 9 + v-V9 = K clt A9 + / K ett = k + k u (2.7) 

where m is a constant depending on the flow properties that can be explicitly calculated 
within the Kraichnan model. Eq. I|2.7|l is just the long-known constant eddy-diffusivity 
closure. 

At the next-leading order we find 

dtd + v-V6 = KZpV a Vpd + f , (2.8) 

with K^g(x,t) — SapK"" — al 2 e a p, and e a p = 1/2 (V a vp + Vpv a ). The filter dependent 
factor a can be determined analytically within the Kraichnan model. Equation i|2.8|) is the 
passive scalar analog (see, e.g. , Kang & Meneveau 2001) of the so-called "mixed-model" 
(nonlinear closure plus scalar eddy-viscosity) used in Navier-Stokes turbulence (see, e.g. , 
Borue & Orszag 1998). Such mixed model was also invoked by Kang & Meneveau (2001) 
to reproduce the correct amount of anisotropy in heated turbulent jets. 
The nonlinear closure can be derived also starting from a Taylor expansion in the spirit of 
Leonard (1974) performed on a modified Leonard term (see, e.g. , Horiuti 1997). A purely 
dissipative effective viscosity model is usually added because the sole nonlinear model 
does not suffice, its dissipation being far too low. Here, the constant eddy-diffusivity plus 
the nonlinear model follow from first-principles without additional requirements. The 
turbulent eddy-diffusivity k t depends on the coarse-grained velocity field: for this reason 
we will dub this model dynamical eddy-diffusivity closure. 

In order to compare the performances of the various models we will explicitly compute 

the structure function (r) = ([9(x + r, t) — 9(x, t)] 2 ) of the coarse-grained scalar field 
according to Eq. 12.7fl or to Eq. (|2.8|) and compare them with the exact value obtained 
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Figure 1. The coarse-grained structure functions obtained from the constant eddy-diffusivity 
closure, and the dynamical eddy-diffusivity model, normalized with respect to the exact filtered 
structure function 13.41 . 



from the solution of Eq. jUJ upon filtering. The result is shown in Fig.^ the dynamical 
eddy-diffusivity closure gives a structure function that is almost indistinguishable from 
the exact result already at scales r w 2 I whereas the constant eddy-diffusivity is not very 
effective in the range r < 10 1. 



3. A systematic approach to LES closure 

3.1. The Kraichnan model 

In order to carry out our analysis, we need to specialize Eq. (|2.1|) to a class of random 
velocity fields. To be more specific, the velocity field is assumed here to be Gaussian, of 
zero mean, statistically stationary, homogeneous and isotropic, (^-correlated in time and 
with inertial-range power-law behavior. Its statistics is fully determined by the correlation 
function 

{[v a { Xl ,t) - v a {x 2 ,t)] Mzi,0) - v p (x3,0)]) = 2DjJ(x x - x 2 )S(t), (3.1) 

where D^)(r) = D a r^[(d + S, — l)5 a p — £r a rp/r 2 ], r = \r\ = \x 2 —x±\ and d is the space 
dimension. The assumption of 5-correlation in time is of course far from the reality, but 
it has the remarkable feature of leading to closed equations for equal-time correlation 
functions C„ = {Q(x\,t) . . .0(x n ,t)) of any order n (see, e.g. , Falkovich et al. 2001). 
The parameter £ governs the roughness of the velocity field, whose Holder exponent is 
£/2. Due to the white-in-time character of the flow, the Kolmogorov value is £ = 4/3. 
A convenient choice for the forcing term is to take / random, Gaussian, statistically 
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homogeneous and isotropic, white in time, of zero mean and with correlation function 

(f( Xl ,t)f(x 2 ,0)) = F(r/L) S(t) (3.2) 

with F (r/L) decreasing rapidly for r L. Since l<CLwe have / ~ /. 

In this framework the expression for the second order correlation function C 2 can be 
derived analytically (Falkovich et al. 2001). In the inertial range rj <C r <C L the exact 
second-order structure function of the scalar field s!f \r) = ([9(xi, t) — 9(x 2 , t)] 2 ) reads: 

where F(0), defined in (|3.2() . is the average injection rate of scalar variance. The exponent 
(2 coincides with the predictions based on dimensional arguments, and is (2 = 2/3 for 
£ = 4/3, according to the Kolmogorov-Obukhov-Corrsin scaling. 

From now on, we will confine ourselves to space dimension d — 3 and specialize Gi to 
a top-hat filter, that is Gi(r) = 3/ (47rZ 3 ) if r < I and otherwise. 

To provide a benchmark for the various closures, we first evaluate the exact value of 
the coarse-grained structure function = ([9{x±, t) — 9(x2, t)] 2 ). A double integration 
on H3.3fl yields 



[^(r)U = Sr« 



(2-0(3-0 fl\ +0 (i 



r I \r 



(3.4) 



Clearly, as the separation r increases and gets far apart from the filter scale I the unfiltered 
result is recovered. Equation l|3.4[) represents therefore the best result that can be achieved 
by means of a closure. 

3.2. Exact statistical equations for unfiltered and filtered fields 

As a first step, we derive the exact equations for the two-point correlation function of 
the filtered and of the unfiltered field. It is more convenient to start the analysis from 
Eq. I|2.1[l by substituting v = v + v' and 9 = 9 + 6' in the advective term v ■ V 9. Equation 
H2.1(l takes then the form: 

d t 9(x 2 ,t) + v(x 2 , t) ■ V 9(x 2 ,t) = K A9(x 2 ,t) + f(x 2 , t) - S(x 2 ,t) , (3.5) 
from which we can immediately derive the equation for C^\r,t) = (9(xi,t)9( X 2,t)}: 

d t C { 2 6 \r,t) + 2{9{x u t)v ■ V9(x 2 ,t)) - 2/coACf \r,t) = F(r) - 2(6(x 1 )S(x 2 )) . (3.6) 

A double convolution of the above equation with the filter Gi yields the exact equation 
for the correlation of the filtered field: 

d t C { 2 6 \r,t) + 2{9(xx,t)v ■ V9(x 2 ,t)) - 2n a AC { 2 \r,t) = F(r) - 2(9(x 1 )S(x 2 )) . (3.7) 

This is the starting point for our systematic procedure to build closure approximations. 
It contains two terms, the second one in the left hand side (l.h.s.) and the last one in 
the right hand side (r.h.s.), which are not expressed as functions of large-scale fields 
only. In the following, we will find approximate closed expressions for the unclosed terms 
perturbatively in l/r. 

3.3. Constant eddy-diffusivity closure 

In order to calculate the second term in the l.h.s. of Eq. 1)3.711 . let us start from (9(xi)S(x2)). 
Its expression can be easily obtained exploiting, e.g. , the Furutsu-Novikov's functional 
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Gaussian integration (see, e.g. , Frisch 1995), which holds for Gaussian velocities and 
forcings. It reads: 

(0( Xl )S(x 2 )) = (4 + 0(6 + {-) --3Q-f [-) +0[~ r ) ■ (3.8) 

Our first claim is that, at order (l/r) 2 , with an error of order (l/r) 2+ ^ we have: 

(6( Xl )S(x 2 )) = (6( Xl )S(x 2 )) = (6(x 1 )S(x 2 )) (3.9) 

and 

(0(xi,t)v ■ V9(x 2 ,t)) = (0(xi,t)v ■ V9(x 2 ,t)). (3.10) 
The derivation of (|3.9|l and l|3.10(l is postponed to the Appendix EI Before proceeding 
further, some comments on H3.10J1 are in order. In plain words, (|3.10|l tells us that the 
Leonard-type term does not contribute, at order 0((l/r) 2 ), to the equation for the second- 
order coarse-grained scalar correlation function. Since our closures are derived from the 
latter equation, it follows that the Leonard-type term will not contribute to small-scale 
parameterizations. This fact is not a consequence of the Kraichnan advection model but 
rather seems to hold for general advection models. Indeed, a standard expansion in the 
spirit of Leonard (1974) (see also, e.g. , Horiuti 1997) performed on (8(xi)C(x 2 )) with 
C given by (|2.5[) . yields at the lowest-order in the filter width the expression 

{9{ Xl )C{x 2 )) ~ l 2 A(9(x x ,t)v ■ V9(x 2l t)). (3.11) 

The latter expression is trivially zero, since (9(x\,t)v ■ W9(x 2 ,t)} is the flux of scalar 
variance which is independent of r = \x 2 — xx\ provided that r falls in the inertial range 
of scales. 

For standard closure models based on single-point quantities, the contribution from the 
Leonard stress in the parameterizations is, generally speaking, nonzero. 

Let us now focus on consequences of relations (13.91) and l|3.10|l on Eq. I|3.7|l evaluated 
at order (l/r)^. At this order, from (|3.8I) and 13.91) we have: 

- 2(6(x 1 )S(x 2 )) = 2/i 1 AC : f ) (r,i). (3.12) 

The above expression follows by comparing (|3.8|) at the order (l/rp with the contri- 
bution coming from the diffusive term, 2k AC^ 0) (r, t) = -F(0)k (3 - £)/(3A)^). We 
immediately realize that the term of order (l/r)* in l|3.8|l corresponds to an effective 
diffusive term with an eddy diffusivity Ki = Z 5 [2«24D ]/[(£ + 4)(£ + 6)]. This gives rise 
to an effective dissipative scale comparable to I: 



(3.13) 



where 770 = (ko/Dq) 1 ^ is the molecular dissipation scale. 

At order (l/r)^, Eq. I|3.7|l is thus closed in the large scales fields and, moreover, due to 
H3.1L)fl it has the same structure of the equation for but with an effective diffusivity, 
Ko + K\. Eq. 12.411 thus takes the form: 

d t 9 + v ■ V 6 = k o11 A9 + f K cff = K + «i- (3.14) 

Starting from (|2.7|) one can deduce the equation for the correlation function 
exploiting the Furutsu-Novikov's functional Gaussian integration. When doing that, the 

next step will consist in comparing the resulting expression for with l|3.4|l . After 
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some simple but quite lengthy algebra, the calculation leads to: 

a(§) _ MO) 

By comparison with l|3.4(l . the above expression permits to quantify the error, which 
occurs at order (l/r) 2 , on the second order structure function due to the closure (|3.14|) . 
The degree of accuracy of the constant eddy-diffusivity description can be perceived by 
looking at Fig. 2] obtained for £ = 4/3, corresponding to the Kolmogorov scaling law for 
the advecting velocity field. 

3.4. Dynamical eddy-diffusivity closure 

Our aim is now to improve the eddy-diffusivity closure, exact at order (l/r)^, by intro- 
ducing a new closure which is accurate up to order (l/r) 2 . It is not difficult, although 
quite lengthy, to verify that the large-scale equation has the form (|2.8I) : 

d t e + i-Ve = Kl V a Vpe + f K^{x,t) = 8 af3 K cii - al 2 e a p , (3.16) 

with e aP = 1/2 {V a vp + W v a ). 

The filter-dependent factor a = J d 3 rGi(r)r 2 /(3l 2 ). For the top-hat filter one immedi- 
ately obtains a — 1/5. 

The equation for (9 2 ), deduced from l|3.16|l . at the statistically stationary state, reads 
(remember that / = /): 

F(0) = 2k^((V6) 2 ) - al 2 {e af} V JV p 6) (3.17) 

which states the energy balance between production (controlled by F (0)) and dissipation. 
The fact that the first term on the r.h.s. of Ij3.17jl gives a dissipative contribution is 
evident. This is actually the case also for the second term. To show that, let us start 
from simple physical considerations. Transforming to the principal coordinates, x' , of 
e Q /3, the term — al 2 e a pSI aOV p6 becomes: 

- al 2 ( Cl \7 v 6\7 v § + c 2 V 2 '0V 2 '0 + c 3 V 3 '0V 3 '0). (3.18) 

Because of incompressibility one has c\ + c 2 + c 3 = 0, and stretching in a given direction 
is always accompanied by compression along, at least, another direction. By virtue of 
the fact that strong scalar gradients are expected to be aligned along the direction of 
maximum compression (corresponding to negative Cj), it then follows that, on average, 
H3.18fl is expected to be positive. This is the mechanism which leads to the well-known 
ramp- and- cliff structure observed in passive scalar turbulence both for Navier-Stokes ve- 
locity fields and for the Kraichnan ensemble (see, e.g., Celani et al. 2001). 
The above arguments can be substantiated within the Kraichnan model. Indeed, exploit- 
ing the Furutsu-Novikov's functional Gaussian integration by parts, one obtains: 

- al 2 {e aP VjV e) = a ^((V^-) 2 )((V0) 2 }. (3.19) 

However, it is clearly possible to observe, locally in space and time, positive values of 
e a /3V Q 0V^0, or, in other words, backscattering events responsible for negative contri- 
butions to the scalar energy flux (see, e.g. , Borue & Orszag 1998 for discussions on 
backscattering events in hydrodynamics turbulence). 

Focusing now on the unfiltered field, the balance equation has the well-known form: 

F(0) = 2k ((V6>) 2 ). (3.20) 



(2-0(3 + (I 



(3.15) 
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Equating the l.h.s.'s of (|3.17|) and Ij3. 20|) . and recalling the inequality n° a + l 2 /Q0((V iVj) 2 ) 
^> kq, one concludes that the gradients of the large scale scalar field are smaller than 
the gradients of the unfiltered field. This is consistent with an effective dissipative scale 
comparable to I and thus much larger than rj by definition. 

To prove Eq. (|3.16() . let us start by rewriting Eq. I|3.7[) with (0(xi)S(x2)) expressed in 
terms of <|3.8[) . Using (|3.9|) we obtain: 

2 



d t C i 2 9) (r,t) + 2(9(x 1 ,t)v- V9(x 2 ,t)) 



2 K " a AC { 2 \r,t) -F(r) = 



F(0K 2 (I 
15 



(3.21) 



where the contribution of order (l/rp in l|3.8l) has been incorporated in the eddy- 
diffusivity term. From l|2.8[) we immediately obtain the equation for C. 



(«). 



d t cP(r,t) +2(6( Xl ,t)v-V6(x 2 ,t)) 

- 2al 2 (d(x 1 )e af3 V a Vf 3 9(x 2 )} 



<9) ( 

0[{l/r) 2 ^] 



2 K Bti AC ( 2 >(r,t) 



F(r) 



(3.22) 



with a = 1/5. We finally have to show that the r.h.s. of (|3.21() and the r.h.s. of l|3.22|l 
coincide up to order (l/r) 2 . In order to evaluate the r.h.s. of l|3.22|l we have to exploit 
again the Furutsu-Novikov's functional Gaussian integration by parts. One thus needs to 
compute the functional derivative S9(x,t)/5v(x",t") which can be easily obtained from 
Accounting for the delta-correlation in time and utilizing the expansions 



D { :hr) = D%(r){l + 



Sjp (r)=Sr(r){l + 



(6), 



(3.23) 



one ends up exactly with the r.h.s. of l|3.21|l . Exploiting once more the Furutsu-Novikov's 
functional Gaussian integration by parts, it is not difficult (although quite lengthy) to 

verify that the expression l|3.4|l for S 2 8 ^ is obtained from 1)2.8(1 . The remaining error is 
pushed at order (l/r) 4 . 



4. Conclusions and perspectives 

Summarizing, a systematic procedure to derive closed dynamical equations for a coarse- 
grained passive scalar field in the statistical steady-state has been obtained in the frame- 
work of the Kraichnan advection model. 

The question that naturally arises, is whether those results are relevant to realistic 
advection models. The answer is given by the outcome of the procedure itself. We recover 
from first principles two well-known closures that are commonly used in applications: the 
constant eddy-diffusivity parameterization of small-scales, and the passive scalar version 
of the nonlinear eddy-viscosity closure used in hydrodynamic turbulence. Of course, the 
value of the effective diffusivity k d " and of the numerical parameter a that appear in 
these closures can be analytically computed only in the Kraichnan model. However, we 
believe that the form of the parameterization can be exported without modifications to 
real situations as well. Clearly, in this case the free parameters (e.g. K off and a) have to 
be determined a posteriori by some empirical procedure. The validity of this approach 
can be checked by direct numerical simulations. 

Let us conclude by mentioning a possible generalization of our work. Our analysis has 
been carried for the second-order correlation function of the scalar field. There are two 
reasons for this choice. First, the second order correlation function is the Fourier trans- 
form of the spectrum of scalar variance, a widely used statistical indicator to characterize 
most of the statistical properties of scalar turbulence. Second, for the Kraichnan model 
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only the second-order correlation function has a simple, closed analytical expression. For 
higher-order correlation functions only perturbative expressions (for example in the limit 
of small £) are available (see Falkovich et al. 2001). However, should have we focused on a 
higher-order correlation function, how would our results change ? Although the analysis 
appears much more cumbersome than the one presented here, the procedure described 
in §|3 can be completed as well: it is still possible to obtain a closed equation for the 
coarse-grained correlation function at any order in l/r, from which one can identify the 
corresponding dynamical equations for the large-scale scalar field. The question is: will 
the latter dynamical equation have the same structure of the coarse-grained scalar equa- 
tion derived from the second-order correlation ? And if this is the case, will the coefficients 
be the same ? Even if the functional form of the closure is preserved, a modification of 
the effective coefficients would mean that strong small-scale fluctuations - associated to 
higher-order correlation functions - must be described by parameters different from the 
ones used for less intense fluctuations. That would question the applicability of closure 
models to the description of the statistics of turbulent fields as temperature or concentra- 
tion, which are characterized by a wide range of fluctuation intensities. This challenging 
issue is left for future research. 

This work has been supported by Cofin 2001, prot. 2001023848 and by HPRN-CT- 
2002-00300 "Fluid Mechanical Stirring and Mixing" . We acknowledge useful discussions 
with Alessandra Lanotte. 



Appendix A. Proof of Eq. (f!T5|) and Eq. (|STTTJ|> 

Let us write, from 1)3.8)1 . (9(xi)S(x 2 )) = A(l/r)^ + B(l/r) 2 . By direct calculation it is 
easily checked that, for any filter Gi (normalized and isotropic), the following relations 
hold up to the second order in l/r. 



(9( Xl )S(x 2 )) 



A 



d 3 sG l (s)(9(x 1 +s)S(x 2 )} 



I 



B 



I 



O 



«+2 



Al- 

. r 



dss z Gi(s) I d(costf) 



2tt 



l-£costf~ + - 



J dss 2 G l (s) J d(cos?9) J df l-2cos??~ + o(^) 

?+2 



Al- 

r 



Ol- 

r 



«+2 



B ' r 



Ol- 

. r 



= (e(x 1 )S(x 2 ))+0[-\ 



iV +2 



Similarly, convolving also over x 2 , one finds: 



(6(xi)S(x 2 )) = (0(x 1 )S(x 2 )) +Ol-\ = (9( Xl )S(x 2 )) +Ol- 



?+2 



«+2 



(Al) 



To prove H3.10J) . let us now consider 1)3.7(1 and the equation obtained from 1)3. 6J) with 
the sole convolution over x\. Because of stationarity all time derivatives vanish and for 
r in the inertial range the two-point terms proportional to molecular diffusivity kq are 
negligible. We can thus rewrite the above two equations as: 



2{6(x x )v ■ V6{x 2 )) = F(r) - 2(0(x 1 )S(x 2 )) 



(A 2) 
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20(x 1 )v ■ V0(x 2 )) = F(r) - 2(0(sc 1 )S(se a )) • (A 3) 

Subtracting (1A~3|) from JXjJl, and using [|3~§j) . we conclude that (0£) = 0[(Z/r)*+ 2 ]. This 
proves (|3.10|) . 
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